Dark energy and the mass of the Local Group 



A. D. Chernin ^'^ 
Received / Accepted 



P. Teerikorpi ^ M. J. Valtonen ^ G. G. Byrd ^ 
L. M. Domozhilova ^ 



V. P. Dolgachev ^ and 



ABSTRACT 

Dark energy must be taken into account to estimate more reliably the amount of dark matter 
and how it is distributed in the local universe. For systems several Mpc across like the Local 
Group, we introduce three self-consistent independent mass estimators. These account for the 
antigravity effect of dark energy treated as Einstein's cosmological constant A. The first is a 
modified Kahn-Woltjer model which gives a value of the Local Group mass via the particular 
motions of the two largest members, the Milky Way and M31. Inclusion of dark energy in 
this model increases the minimum mass estimate by a factor of three compared to the classical 
estimate. The increase is less but still significant for different ways of using the timing argument. 
The second estimator is a modified virial theorem which also demonstrates how dark energy 
can "hide" from detection a part of the gravitating mass of the system. The third is a new 
zero-gravity method which gives an upper limit to the group mass which we calculate with high 
precision HST observations. In combination, the estimators lead to a robust and rather narrow 
range for a group's mass, M. For the Local Group, 3.2 < M < 3.7 x lO^^M©. Our result agrees 
well with the Millennium Simulation based on the ACDM cosmology. 

Subject headings: galaxies; Local Group; cosmology: dark matter, dark energy 



1. Introduction 

It has long been known that galaxy groups and 
clusters would fly apart unless they were held to- 
gether by the gravitational pull of much more 
material than what is actually seen. The rapid 
progress in dark matter studies stems from opti- 
cal surveys of large areas and high redshifts, CMB 
fluctuation measurements, sharp X-ray images, 
gravitational lensing measurements, etc. Concur- 
rently, conclusive observational indications have 
accumulated that dark matter and baryons con- 
tribute not more than 30% to the mass/energy 
content of the universe. The rest 70% is composed 
of a much more mysterious dark energy (Riess et 
al. 1998, Perlmutter et al. 1999). Dark energy 
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does not cluster, and its simplest mathematical 
formulation goes back to 1917 when Einstein in- 
troduced the cosmological constant A. This in- 
terpretation is assumed in the currently standard 
ACDM cosmology. A positive A corresponds to a 
positive constant energy density in all entire space. 
Such a medium with negative pressure (pi, > 0, 
Py < 0) is characterized by the equation of state 
Py ~ — py (c = 1), like that of vacuum. According 
to Einstein's equations, gravity depends on pres- 
sure as well as density: the effective gravitating 
density pcs = p + ip is negative for a vacuum 
{— — 2p^), leading to a repulsion, or antigravity. 
The dark energy antigravity revealed itself first in 
the accelerating cosmological expansion. 

In this paper, we follow the ACDM cosmology 
and assume that the dark energy (DE) is a perfect 
fluid with a constant density every where and in 
any reference frame. From most recent CMB stud- 
ies (Spergel et al. 2006), p„ = (0.72±0.03) x lO'^s 
g cm~'^. We focus on relatively small spatial 
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scales, on the seale of groups of galaxies. We show 
that galaxies and their systems "lose" a part of 
their gravity due to the antigravity of the dark en- 
ergy in their volumes. The dynamical methods of 
the mass estimation of dark matter (and baryons) 
should take into account this "lost-gravity" effect 
- otherwise the mass would systematically be un- 
derestimated. 

2. Modified Kahn-Woltjer (MKW) esti- 
mator 

Fifty years ago, Kahn & Woltjer (1959, KW) 
used simple linear two body dynamics to describe 
the relative motion of the Milky Way and AI31 
galaxies. The motion of the galaxies was described 
(in the reference frame of the binary's center of 
mass) by the equation of motion 

D{t) = -GM/D^, (1) 

where D is the distance between the galaxies and 
M is the total mass mi + of the binary. In the 
well-known first integral 

^-b{tf = GM/D + Eq, (2) 

the constant total mechanical energy of the binary 
(per unit reduced mass) is negative for a gravita- 
tionally bound system, £^0 < 0- This fact and 
Eq.2 with the observed values D = 0.7 Mpc and 
D = —120 km/s lead to an absolute lower limit 
for the estimated binary mass: M > 1 x lO^^M©. 

The "timing argument" KW calculation inte- 
grates in reverse the equation of motion (1) from 
the above present-day separation and radial ve- 
locity requiring that the initial expansion of the 
two galaxies from one another started a time in 
the past equal to the age of the universe (13.7 bil- 
lion years). Thus we calculate a mass for the pair 
of M ~ 4.0 X 10^2 Af0. See Binney & Tremaine 
(1987) for an analytic discussion (for the 10 to 20 
billion year age range suspected at that time). 

With a minimal modification of the original 
KW method, including the dark energy back- 
ground pv, the equation of motion and its first 
integral become: 

D{t) = -GM/D"^ + (3) 
^-b(tf = GM/D + G^p.D^ + Eo. (4) 



Now the total energy for a bound system em- 
bedded in the DE background is 

Eo < -^GMV3(|[p,)i/3, (5) 

as may be easily seen from the gravitational po- 
tential U = -GM/D - ^Gp^D^ in the right side 
of Eq.4. With the same data as above for current 
separation and the relative velocity of the galaxies, 
Eqs. 4,5 lead to a new absolute lower mass limit: 
M > 3.2 X IO^^Mq, or 3 times the absolute lower 
limit of the original statement of the problem. 

Again the MKW equations can be integrated 
back to the origin of the universe including the 
dark energy value given earlier. The same data 
lead to a larger mass: M = 4.5 x 1O^^M0. 

In fact, the KW model can hardly be extrapo- 
lated to the origin of the universe, while it seems 
realistic for the stage of the binary collapse. In 
Fig. 1 we show the calculated mass as a function 
of the time back to the maximal separation of the 
Milky Way and M31; the lower and upper curves 
show the results without and with dark energy, re- 
spectively. In the ACDM cosmology gravitational 
instability is terminated in the linear regime about 
7 Gyr ago with the start of the DE dominated 
epoch (e.g. Chernin et al. 2003). Taking this 
to be the minimal collapse time yields the maxi- 
mal mass of the binary, and this is 3.3 x 
without dark energy and 4.1 x lO^^M© with dark 
energy. 

These results demonstrate clearly the effect 
of the lost gravity in gravitationally bound sys- 
tems embedded in the dark energy background. 
The relative motion of the bodies is controlled by 
the gravity (of dark matter and baryons) which 
is partly counterbalanced by the DE antigrav- 
ity. Consequently, the mass estimate must be cor- 
rected: the mass of the Local Group given by the 
modified Kahn-Woltjer (MKW) estimator is sig- 
nificantly larger than that obtained without in- 
clusion of dark energy. Here we have taken ad- 
vantage of our position within the Local Group, 
the well known distance of M31 and the measured 
approach speed. We thus should modify conven- 
tional ways to estimate the masses of groups of 
galaxies in general. 
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3. Modified virial theorem (MVT) 

Conventional virial mass estimators use the re- 
lation between the mean total kinetic < K > and 
potential < U > energies of a quasi-stationary 
gravitationally bound many-body system: < 
K >= 1/2 1 < U > \. The presence of dark 
energy modifies this relation. The total potential 
energy includes now not only the sum Ui oi the 
mutual potential energies of its member particles, 
but also the sum U2 of the potential energy of the 
same particles in the force field of dark energy: 



1 ^ Grrnmj 



4iTpy 2 

U2 = ^2^mi^^ 



(6) 

Here rj is the radius-vector of a particle in the 

frame of the system's barycenter; the summation 
in Ui is over all particle pairs (i 7^ j). The major 
contribution to the sum is from dark matter par- 
ticles whatever their individual masses may be. 

A hint to the structure of U2 may be seen, e.g., 
from the second item in the right side of Eq.4; the 
summation in U2 is over all the particles. Dark 
energy comes to the virial theorem via an extra 
contribution to the potential energy of the system. 

A link to the kinetic energy is provided by the 
equation of motion of an individual particle: 



dU dUi 
rriiTi = -— = -—— 



6 n 



(7) 



Averaging over time and using the Euler theorem 
on homogeneous functions (applied separately to 
the two functions that come from the two terms 
in the right side of Eq.7), we find 



< K >= -- < Ui > + < U2 > . 



(8) 



This is the new virial relation adapted to the uni- 
verse with the dark energy background. 

Eq.8 may be rewritten in terms of the total 
mass M, a characteristic velocity V and charac- 
teristic sizes Ri,R2: 



M 



V^Ri Stt -3 



(9) 



or in convenient units: 

^=2.3xl0«(-^, 
Mq ^km/s^ 



)2(JL) + 0.9x1012(^) 



R2 



'Mpc 



Mpc' 
(10) 



This new mass estimator - the modified virial the- 
orem (MVT) estimator - includes an additional 
positive term which is equal to the absolute value 
of the effective (anti)gravitating mass of dark en- 
ergy contained in the spherical volume of the ra- 
dius i?2: Afe// = —^PvR2- This term gives a 
quantitative measure of the lost-gravity effect. 

One can get a physical sense of Eqs.9,10 by us- 
ing representative values of the characteristic size 
.Ri, R2 and velocity V. In the simplest example of 
only one body orbiting a gravitating mass M, the 
size .Ri = R2 is the radius of the orbit, the veloc- 
ity V is the orbital velocity. In this case, the total 
effective gravitating mass within the orbit is the 
gravitating mass M plus the antigravitating (neg- 
ative) mass Meff. A similar identification of the 
quantities is obvious as well when the dark mass 
of a galaxy is derived from its rotation curve. In 
both cases, the estimated mass is larger than that 
in conventional estimates: the additional mass is 

\Meff\. 

Generally, the characteristic sizes and veloci- 
ties for groups and clusters need more sophisti- 
cated analysis - even in the absence of dark en- 
ergy (e.g. Peebles 1971). Here we merely point 
out the expected significance of this effect in small 
systems like the Local Group. The first (conven- 
tional) term in Eq.9 is estimated for the Local 
Group as 2.3 x IO'^'^Mq (van den Bergh 1999) or 
1.9 X lO^^M© (Karachentsev 2007, Karachentsev 
et al. 2009). Then the second term in Eq.9 in- 
creases the total mass by 30 to 50%, if the total 
size R2 is about 1 Mpc. Such an effect is compat- 
ible with the other methods here discussed. 

Assuming _Ri = R2, one may see that the rela- 
tive contribution of the dark energy scales as the 
crossing time squared {R/V)"^. If so, the lost- 
gravity effect is 10-30 times larger in a group like 
the Local Group than in a rich cluster like the 
Coma cluster. 

4. Zero-gravity surface (ZGS) estimator 

As it was mentioned elsewhere (Chernin 2001, 
2008; Byrd, Chernin & Valtonen 2007; Teeriko- 
rpi et al. 2008), an isolated gravitationally bound 
galaxy group can exist on the dark energy back- 
ground, only if gravity dominates over antigravity 
in the whole volume of the system. In a simple 
model, a group may be represented by a spheri- 
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cal mass M of dark matter and baryons, embed- 
ded in the uniform dark energy background. Out 
of the mass at a distance R from the mass cen- 
ter, the gravity force is given by Newton's inverse 
square law in the reference frame related to the 
group barycenter. The antigravity force produced 
by the DE density is given by Einstein's linear 
law (e.g. Chernin 2001, 2008, Chernin et al. 2006) 
in the same reference frame: 

F^ = -GM/R^, Fe = +yGpvR. (11) 

Gravity and antigravity are exactly balanced at 
the zero-gravity surface of the radius R = Ry 
(Chernin et al. 2000; Chernin 2001; Baryshev et 
al. 2001, Dolgachev et al. 2003): 

Rv = i^y^\ Rv>Ro, (12) 

where Rq is the radius of the group. Antigravity 
dominates at i? > i?„ , and gravity is stronger than 
antigravity a,t R < Ry. The lost-gravity effect 
is one hundred percent value on the zero-gravity 
surface: the total effective gravitating mass of the 
non- vacuum matter and vacuum energy contained 
within the surface is zero on the surface. No fi- 
nite bound orbits are possible on the surface and 
further at distances R> R^. 

Since the zero-gravity radius Ry depends on the 
gravitating mass of the group alone (for a fixed 
p„), this mass may be found, if the zero-gravity 
radius is known from observations: 

M = y ~ 0.9 X 10i2[i?^/(Mpc)]3MQ. (13) 

The position of the zero-gravity surface can 
indeed be detected. The most complete high- 
precision data come from systematic observa- 
tions of galaxies up to the distance 3 Mpc per- 
formed recently with the Hubble Space Tele- 
scope (Karachentsev 2005, 2007; Karachentsev 
et al. 2002, 2003, 2009). About 60 galaxies with 
good distances (8-10% accuracy) and velocities 
(5-10 km/s accuracy) are observed in the volume. 
Nearly half of them join the Milky Way and M31 
galaxies forming the Local Group. In the velocity- 
distance diagram they occupy the distance range 
of about 1 Mpc; their radial velocities (relative 
to the group barycenter) are in the interval from 
-150 to -1-150 km/s. The other half of the ob- 
served galaxies (all of them are dwarfs) move away 



from the group and have only positive velocities 
in the diagram. It is important that their min- 
imal distances are 1.1-1.6 Mpc from the group 
center (Karachentsev 2005. Chernin et al. 2007, 
Karachentsev et al. 2009). The flow of receding 
galaxies is rather regular: it follows mainly the 
linear velocity-distance relation, and its velocity 
dispersion is only about 25-30 km/s. 

This radical difference in the phase-space struc- 
ture of the group and the outflow around it is a 
most remarkable feature. We argue (Chernin et al. 
2000, 2007, Chernin 2001, Teerikorpi et al. 2008) 
that the physics behind this feature might be due 
to the interplay between the gravity of the group 
and antigravity of the dark energy background, so 
that the group size Rq is less than R^ and the 
flow starts at the distances R > Ry. Therefore 
the zero-gravity surface is located somewhere in 
the gap between the group and the outflow, i.e. in 
the distance interval from 1.1 to 1.6 Mpc. 

The double inequality 1.1 < Ry < 1.6 Mpc and 
Eq.l3 lead to the lower and upper limits for the 
non- vacuum mass of the group: 

1.2 < M < 3.7 X lO^^M©. (14) 

This result is compatible with the estimations 
obtained in Sees. 2, 3. The ZGS-estimator uses the 
overall structure of the group as a whole, while 
the MKW- and MVT-estimators deal with the 
internal dynamics of the group. The new ZGS- 
estimator is possible only due to the existence of 
dark energy in the universe. 

5. Conclusions and discussion 

The local volume proves to be a prospective 
arena for the study of dark matter and dark en- 
ergy (e.g. Chernin 2001, 2008; Byrd, Chernin & 
Vahonen 2007; Niemi et al. 2007; Teerikorpi et 
al. 2008). We have demonstrated above that, sur- 
prisingly enough, by taking the dark energy back- 
ground into account enables one to learn more on 
the mass of dark matter in galaxies and their sys- 
tems. 

To summarize the results: 

1. The minimally modified Kahn-Woltier 
method leads to a new absolute minimum of the 
Local Group mass: M > 3.2 x IO^^Mq, and it 
is 3 times larger than the similar mass minimum 
obtained via the original KW method. The back 



4 



integration leads to a mass of the M31, Milky Way 
pair of 4.5 x lO^^M©, significantly greater than the 
original method. 

2. The modified virial theorem has the form: < 
K >= -5 < C/i > + < C/2 >• It accounts for the 
potential energy U2 which is due to the antigravity 
force of dark energy. The theorem gives a clear 
measure of the lost-gravity effect. 

3. The Local Group is located in the area of 
its self-gravity domination which makes it grav- 
itationally bound, while the outflow around the 
system develops in the area where no bound or- 
bits are possible. The zero-gravity surface sepa- 
rates these two areas, and its location identified 
with the observed velocity-distance diagram leads 
to the absolute upper limit of the group mass: 
M <3.7 X IQI^Mq. 

Some additional remarks: 

1. The linear binary model for the relative mo- 
tion of the Milky Way and M31 galaxies (Sec.2) 
is an idealization. The real dynamics of the Lo- 
cal Group must be more complex (see the study 
by Valtonen et al. 1993). In particular, a trans- 
verse velocity might increase the mass value, but 
not more than by 10-15%. The earlier dynami- 
cal history of the binary (not even mentioning its 
formation) can hardly be described by this model. 
This process needs more studies, and cosmological 
N-body simulations are a real tool to clarify the 
matter. The most recent Millennium Simulation 
shows a mass range (2 — 5) x IO^^Mq for galaxy 
groups similar to the Local Group (Li and White 
2008); it seems instructive that the interval found 
in Sees. 2,4 is in a good agreement with the simu- 
lation result. 

2. A study of the lost-gravity effect in groups 
in general is complicated by the fact that in ex- 
isting group catalogues, the virial masses tend to 
be overestimated for other reasons (Niemi et al. 
2007); we will discuss this in detail elsewhere. We 
will also reexamine the large-scale mean dark mat- 
ter and baryonic densities and study the related 
problem of the dark mass deficit in the local vol- 
ume (Karachentsev et al. 2009). 

3. The gravity potential of a galaxy group 
like the Local Group which contains a dominant 
binary of giant galaxies is obviously non-static 
and non-spherically symmetrical - contrary to the 
model of Sec. 4. However non-sphericity and time- 



dependence are not significant when one considers 
relatively large distances from the group barycen- 
ter. Calculations show that the deviations of the 
zero-gravity surface from sphericity now and dur- 
ing the binary collapse are not larger than about 
10% in radius (Dolgachev et al. 2003; Chernin et 
al. 2004) . The deviations are even less in the area 
of the outflow around the group. 

4. It is worthwhile to compare the zero-gravity 
method of Sec. 4 with the classical method of the 
zero-velocity distance Rq predicted by Lemaitre- 
Tolman solution without or with dark energy 
(e.g. Peirani & de Freitas Pacheco 2006, 2008, 
Karachentsev et al. 2009). In that approach the 
equation of motion for mass shells surrounding the 
central mass M is solved by integrating from a 
moment near the Big Bang. The resulting present 
velocity-distance relation gives the turn-around 
radius Rq, to be compared with the observations. 
The mass turns out to be proportional to H^Rq 
and the coefficient of proportionality is different 
for the solutions with or without the cosmological 
constant, so that for a given Rq the inferred mass 
is again larger when dark energy is included (by 
about 30%; Peirani & de Freitas Pacheco 2006). 
In contrast, our work on the early dynamics of the 
Local Group suggests that there is a finite time 
interval after the Big Bang where the outflow 
of dwarf galaxies originates (e.g. Chernin et al. 
2004, 2007; Byrd, Chernin & Valtonen 2007). In 
this case the present zero- velocity distance will not 
be uniquely predicted, while, on the contrary, the 
zero-gravity distance has its near-constant value 
depending just on M and py. 

Peirani & de Freitas Pacheco (2008) did not 
use directly Rq, but fitted their LT-predictions 
through velocity-distance data points around the 
Local Group and some nearby groups. They did 
not get different mass estimates with or withoiit 
the cosmological constant. A possible reason may 
be the assumption of the simultaneous origin of 
the outflow near the cosmological singularity, and 
their value 2.4 x lO^^M© should be regarded as a 
lower limit. 

To conclude, the mass determination for galax- 
ies and their systems is still a hard observational 
problem - even for the nearest well studied objects. 
The task needs also a reliable theory background. 
The lost-gravity effect is a new significant element 
of the underlying physics. It is included in three 
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different and self-consistent mass estimators intro- 
duced above. Their results are compatible with 
each other, and in combination they give a robust 
and rather narrow range of the Local Group non- 
vacuum mass: 3.2 < M < 3.7 x lO^^Mg. Even 
for the lower limit (from the MKW binary energy 
model), the individual dark halos of the Milky 
Way and M31 galaxies contain together at least 
90% of the total dark matter of the group, and 
the rest might be distributed in the outer common 
halo of the group around its dominated binary. 

We thank A. Cherepashchuk, Yu. Efremov, A. 
Silbergleit and A. Zasov for discussions. AC ap- 
preciates a partial support from RFBR under the 
grant 09-02-00418. 
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Fig. 1. — This diagram shows the difference between the Local Group mass predictions from the classical 
Kahn-Woltjer estimator (KW) and its modified form here introduced (MKW). The x-axis gives the look-back 
time of maximum separation and the y-axis gives the calculated total mass of the Milky Way & M31 binary. 
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